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1. Introduction . Current theories, in which the observed antiproton component is 
attributed strictly to secondary production in high energy inelastic collisions of 
protons with the interstellar medium or the atmosphere, apparently fail to explain 
the relatively high p vertical intensities measured at mountain and balloon 
altitudes. Therefore, a more careful calculation of the theoretical secondary 
intensity spectra is required before more exotic sources for these excess p's can 
be explored. 

In this paper, we have used a one-dimensional diffusion equation (valid if 
®lab - 20° down to sea level) to calculate the expected vertical intensity of p's 
due only to secondary production in the atmosphere; any assumed primary p 
spectrum can also be included. Two adjustable parameters, the inelasticity and 
charge exchange in nucleon-nucleus collisions, were included in our algorithm. In 
order to obtain an independent estimate of their values, we first calculated the 
proton vertical intensities in the atmosphere, adjusting the parameters until our 
curves fit the experimental proton data, and then assumed that these values were 
identical in antinucleon-nucleus collisions. 


2. Results . Our calculations followed a method suggested to us by T. K. Gaisser 
in which the atmosphere was divided into ’’slabs" of equal thickness Ax; slabs of 1 
g/cm 2 were used. In calculating the differential proton intensities, a primary 
proton spectrum of the form 1 j p = 2(E + 2.15)" 2 * 75 cm" 2 -s’ 1 -sr“ 1 -GeV' 1 , where E 

is the proton kinetic energy in GeV, was assumed. Protons and neutrons from 

higher Z nuclei were assumed to have the same spectral shape, and all protons and 
nuclei whose momenta were less than the vertical cutoff rigidity were excluded. 
Then, working from the top of the atmosphere down to sea level, the proton 
intensity of the i + 1 slab was calculated using the equations n(i + 1) = n(i) + 
(dn(E)/dx)Ax and 


dn p (E) 

dx 


^ + («)^B + T 1 dN(E,E Q ) 

.air dx Ae r _p.air dE 
A. . tn _tL A. . 

inel u inel 


[(l-a)n p (E 0 ) + an n (E 0 )]dE 0 , 


( 1 ) 


with a similar equation [without the ionization loss 2 term (dE/dx)(An p /AE)] for the 
neutron intensity, n p . All of the values of n on the right-hand side are the values 
from the slab i above. is the inelastic mean-free path of proton-air nuclei 

collisions, scaled from p- 1 C data. 3 ’ 4 The last term in the equation above adds in 
the protons gained from inelastic collisions of higher energy protons with air 
nuclei, with dN(E,E 0 )/dE defined as the probability of a proton with initial energy 
E 0 possessing energy E after collision. A uniform probability distribution ranging 
from 0 to eE 0 for dN(E,E 0 )/dE, with average elasticity e/2, was assumed, e and 
the probability a of charge exchange were our adjustable parameters. With the 
values a = 0.333 and e = 0.9, our computed proton intensities matched the 
experimental data quite well (see Fig. 1). We then used these values for e and a 
in our antiproton intensity calculations. 
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Fig. 1. The curves show the 

calculated results for vertical proton 
intensities at 710, 747, and 1030 
g/cm 2 depth. The data is a 
compilation by Barber et al., Ref. 14. 


Using the same method, the antiproton intensity due only to secondary 
production was calculated with eguations analogous to Eg. (1), but with the term, 


daJ ir (E,E 0 ) 


air 
a. i 
inel 


n nucleon^0^0 


added on the right side for the production spectrum of p's in proton-air collisions, 
which is taken from a parameterization of Tan and Ng, 5 attenuated with an 
attenuation length of 122 g/cm 2 as depth increases. The form of dE/dx and 
dN(E,Eg)/dE were the same as used in the proton calculations. Xj^ lr , the p mean- 
free path in air for annihilation and inelastic scattering, was estimated by scaling 
a power-law fit to p- 12 C reaction cross sections from data compiled by G. Bruge 8 
and provided to us by J. C. Peng, LAMPE; the result is shown in Fig. 2, curve d, 
along with another estimate, 7 curve e. Three different forms of the 

antiproton mean-free path in air for inelastic, non-annihilation collisions* were 
tried, as shown in Fig. 2. In version c, X^ir ^ W as assumed to be equal to for 

protons. A curve from Szabelski and Wolfendale 7 was used in version b. Fjnally, 
in version a, we attempted a realistic estimate for Af^l by assuming that Xf^^j / 
X§ffnih = a annih^PP^ a inel^PP^ 5ince there is little data available on a inel (pp) [non- 
annihilation], we assumed that it is the same as the pp inelastic cross section. 8 
Below 0.5 GeV, X^i^i depends entirely on quasi-elastic p-nucleon scattering; for 
our realistic estimate shown in curve a, the quasi-elastic scattering was taken to 
be 1/10 as probable as for the p-nucleon case, based upon special Monte Carlo 
runs using ISABEL INC for p- 2 C inelastic scattering at 180 and 400 MeV 9 
arranged for us by P. L. McGanghey at Los Alamos. 
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Fig. 2. Antiproton-air interaction 
lengths (mean-free paths) em- 
ployed in the calculations: (a) 

— air 

Xj ne l derived from our most real- 
istic estimate of aj ne j(p-air); (b) 

air -air 

Ajnel from Ref. 7; (c) Ai ne i derived 

if a inel(P- air ) = a R (P- air )» (d) Xr 11 " 
derived from our fit to a R (p- l 2 C) 

— gjp 

data; (e) ar from Ref. 7. 


The results of our antiproton calculations, along with experimental data for 
the p vertical intensities at mountain altitudes, are shown in Fig. 3. At mountain 
altitude, both curves a and c are consistent with observations by the Arizona 
group, 1 ® although we feel that the comparison with curve a, our most realistic 
estimate, is the more significant one. At balloon altitudes, our secondary p 
intensity estimates are an order of magnitude below the intensities reported by 
Buffington et al . 11 and Golden et al . 12 




Fig. 3. The curves show calculated antiproton intensities produced 
by atmospheric interactions at ^47 g/cm 2 11 g/cm 2 , and 5 g/cm 2 
depth. The data points are taken from Refs. 10 , 11 , and 12. The 
designations (a), (b), and (c) correspond to utilizing curves (a,d), (b,e), 
and (c,d), respectively, from Fig. 2. 
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These results assume that there is no primary p spectrum. We then 
added to version "a" a primary p spectrum normalized by passing through 
the pointy of Golden et al. 12 at 8.5 GeV of the form 

then 
) was 

2.0% greater than the purely secondary p spectrum; with y-2A as 
suggested by Stecker and Wolfendale, 13 the primary p contribution at 
mountain altitude increases to 4.8%. Then y was adjusted until the 
calculated low energy p intensity at 747 g/cm 2 , due to the primaries, 
equaled the difference between the data from Bowen et al. 10 and the 
calculated result for purely secondary p spectrum [Fig. 3(747 g/cm 2 ), curve 
a]. A rather flat spectrum, y=0.95, is required. We also determined the 
most probable energy of the primary antiprotons contributing to the low 
energy p's at 747 g/cm 2 : for y=2.1, it was 30 GeV; for y=1.5, 240 GeV; for 
y=0.9, 3800 GeV. Obviously, more data on the antiproton intensities at high 
altitudes, as well as additional data on p cross sections, are needed before 
making an analysis with fewer approximations. 
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